Topological Defects and Interactions in Nematic Emulsions 
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Inverse nematic emulsions in which surfactant-coated water droplets are dispersed in a nematic 
host fluid have distinctive properties that set them apart from dispersions of two isotropic fluids 
or of nematic droplets in an isotropic fluid. We present a comprehensive theoretical study of the 
distortions produced in the nematic host by the dispersed droplets and of solvent mediated dipolar 
interactions between droplets that lead to their experimentally observed chaining. A single droplet 
in a nematic host acts like a macroscopic hedgehog defect. Global boundary conditions force the 
nucleation of compensating topological defects in the nematic host. Using variational techniques, 
we show that in the lowest energy configuration, a single water droplet draws a single hedgehog 
out of the nematic host to form a tightly bound dipole. Configurations in which the water droplet 
is encircled by a disclination ring have higher energy. The droplet-dipole induces distortions in 
the nematic host that lead to an effective dipole-dipole interaction between droplets and hence to 
chaining. 

PACS numbers: 77.84.Nh, 61.30.Cz, 61.30.Jf 
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I. INTRODUCTION 

Topological defects Q-^] , which are a necessary conse- 
quence of broken continuous symmetry, exist in systems 
as disparate as superfluid helium 3 01 and 4 ||, crys- 
talline solids liquid crystals ]lO|^d| , and quantum- 
Hall fluids ]12| |. They play an important if not deter- 
mining role in such phenomena as response to external 
stresses ||,|| , the nature of phase transitions [||,[l3|,|l4| , or 
the approach to equilibrium after a quench into an or- 
dered phase [ fl5|| ; and they are the primary ingredient in 
such phases of matter as the Abrikosov flux-lattice phase 
of superconductors ||l6| , |l7f or the twist-grain-boundary 
phase of liquid crystals |l8|-p(| . They even arise in cer- 
tain cosmological models||2l]]7 Given the universal na- 
ture of topological defects, it is always interesting to find 
new systems that allow us to increase our understand- 
ing of these defects. In this paper, we will present a 
detailed theoretical investigation of a new class of ne- 
matic emulsions [B2| whose intriguing properties are con- 
trolled by a class of topological defects called hedgehogs. 
These emulsions are either simple inverse emulsions in 
which surfactant-coated water droplets are dispersed in 
an aligned nematic host, or they are multiple emulsions 
in which water droplets are dispersed in larger nematic 
drops that in turn are dispersed in water. 

Liquid crystals are ideal materials for studying topo- 
logical defects. Distortions yielding defects are easily pro- 
duced through control of boundary conditions, surface 
geometries, and external fields. The resulting defects are 
easily imaged optically. The many different liquid crys- 
talline phases (nematic, cholesteric, smectic-A, smectic- 
C, etc.) with different symmetry ground states make 
it possible to study different kinds of defects. Over the 



years, liquid crystals have provided us with detailed and 
visually striking information about topological defects. 

Liquid crystal emulsions in which surfactant-coated 
drops containing a liquid-crystalline material are dis- 
persed in water have been a particularly fruitful medium 
for studying topological defects [^3|- ^6| , ^0| . The liquid- 
crystalline drops are typically from 10/mi to 50/itm in 
diameter and are visible under a microscope. Changes 
in alignment direction, specified by the Frank director 
n, are easily studied under crossed polarizers. The iso- 
lated drops in these emulsions provide an idealized spher- 
ical confining geometry for the liquid crystal. More gen- 
eral distorted or multiply connected random geometries 
p6| such as those produced in polymer-dispersed liquid 
crystals (PDLCs) p7| , p8| , in emulsion films, or in dis- 
persions of agglomerations of silica spheres in a nematic 
host J|9| are of considerable current interest because of 
display technologies based upon changing the light scat- 
tering properties of these systems through modification 
of defect distributions via external fields. 

In this paper, we will study inverse and multiple ne- 
matic emulsions. These emulsions differ from the di- 
rect emulsions described above in that isotropic water 
droplets are dispersed in a nematic host rather than the 
other way around. They are considerably more complex 
than direct emulsions. In direct emulsions, the nematic is 
separated into distinct, nearly spherical drops. Normal or 
homeotropic boundary conditions on the nematic direc- 
tor at a drop's surface will lead to a single point hedgehog 
defect in its interior; tangential boundary conditions will 
lead to a pair of surface defects called boojums ]30| , p5| , ^l[ . 
Though there can be transitions among various director 
configurations as temperature or boundary conditions are 
changed 32 33] , the topological structure of these drops 



1 



is simple. In inverse emulsions, each water drop with 
homeotropic boundary conditions will create a hedgehog 
director configuration in its immediate vicinity. Global 
boundary conditions at the surface of the nematic restrict 
total topological charge. Thus, in order to satisfy global 
constraints, additional defects must be created out of the 
nematic to compensate for or to cancel the topological 
charge created by droplets. The nature and placement 
of these additional defects determine the far-field direc- 
tor distortion produced by a droplet and the nature of 
droplet-droplet interactions. Experiments [^2|,^4| show 
that each water droplet creates a companion point de- 
fect leading to dipole distortions of the director field at 
large distances. This is in contrast to the quadrupolar 
"Saturn-ring" configuration in which a disclination ring 
encircles a dro plet at its equator that has been exten- 
sively studied |35|-|3ql. Our calculations show that the 
experimentally observed dipole configuration is the pre- 
ferred one and that it leads to a dipole-dipole interaction 
between drops that gives rise to the experimentally ob- 
served chaining of droplets. It is interesting to note that 
similar topological dipole configurations appear in two- 
dimensional systems including (1) free standing smectic 
films fl39| where a circular region with an extra layer plays 
the role of the emulsion water droplet and (2) Langmuir 
films jijfl in which a liquid-expanded inclusion in a tilted 
liquid-condensed region plays a similar role. 

The outline of this paper is as follows. In Sec. II, 
we review important elastic and topological properties of 
nematics. In Sec. Ill, we provide an overview of impor- 
tant experimentally observed properties of inverse and 
multiple nematic emulsions. In Sec. IV, we calculate 
the director configurations and energy of a single water 
droplet in a uniform nematic using various variational 
ansatzes. In Sec. V, we introduce a phenomenological 
free energy to describe long-distance director distortions 
and interactions among droplets. Finally, In Sec. VI, we 
summarize our results. 



II. ORDER, ENERGY, AND TOPOLOGICAL 
DEFECTS IN NEMATICS 

A nematic liquid crystal is a uniaxial, homogeneous 
fluid characterized by a unit vector n, called the Frank 
director, specifying the direction of the principal axis 
of a symmetric-traceless-tensor order parameter. The 
ground-state free energy of a nematic is invariant under 
all spatially uniform rotations of n and under all inver- 
sions n — > n. In addition, all physically observable 
quantities are invariant under n — * — n. The ground- 
state manifold or order-parameter space is the unit 
sphere in three dimensions S 2 with opposite points iden- 
tified, i.e., the projective plane RP 2 = S 2 /Z 2 @,|.|To|. 
The topological structure of the ground-state manifold 
determines the types of possible topological defects. As 
we will review below, nematics can have both line defects 
(disclinations) and point defects (hedgehogs). 



A. The Frank Free Energy 

The energy of slowly varying spatial distortions of the 
director n(r) is determined by the Frank free energy 

F 1 = | J d 3 r{ifi(V • n) 2 + K 2 (n • V x n) 2 

+K 3 [nx (V x n)] 2 } (1) 
- J d 3 rA 24 V • [n x (V x n) + n(V • n)], 

where Ki, K2, K3, and A24 are, respectively, the splay, 
twist, bend, and saddle-splay elastic constants. (There 
is also the possibility of another surface term with en- 
ergy if 13V • (nV • n) |il]j42]| , which we will not consider 
in this paper.) The saddle splay term is a pure diver- 
gence; it reduces to integrals over all surfaces, including 
interior surfaces formed, for example, by water droplets. 
For spherical surfaces with normal boundary conditions, 
these integrals are constant and do not vary, for exam- 
ple, when separations between droplets are changed. To 
keep our calculations as simple as possible, we will use 
the one-constant limit of the Frank free energy: 



F = \K J d 3 r[(V • n) 2 + (V x n) 2 ] 



(2) 



- A 24 J dS • [nV • n + n x (V x n)] 

= \K J d 3 rV l n J \7 i n J (3) 

+\{K - 2A 24 ) J dS ■ [nV • n + n x (V x n)]. 

Since surface energies do not play an important role in 
the phenomena to be discussed in the paper, we will set 
the saddle-splay constant A 24 equal to zero unless other- 
wise specified. When K — 2A 24 , the free energy reduces 
to the first line of Eq. (||), which is invariant with respect 
to rigid rotations of any director configuration. [Note: In 
Ref. p2| , calculations were done with K = K24.} 

B. Surface Energies 

Surfaces generally impose a preferred alignment direc- 
tion of the nematic director relative to their local nor- 
mals. The energetics of this alignment are described by 
the Rapini-Papoular phenomenological surface free en- 
ergy pi 



F S = W I dS"sin 2 7 



(4) 



where 7 is the angle between the director and the surface 
normal. Homeotropic or normal alignment is favored by 
W > and tangential alignment by W < 0. The coupling 
constant W varies in the range 10 _4 -1 erg/cm 2 [Q with 
typical values of order 10~ 3 erg/cm 2 pa]. 
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In addition to the above surface alignment energy, in 
emulsions, there is the energy arising from the surface 
tension of the water-surfactant-oil interface. This energy 
is simply the surface tension a times the total surface 
area 



dS. 



(5) 



The surface tension is of order 10erg/cm 2 46 1. 

We can now discuss the relative importance of the 
surface energies and the bulk Frank energy. Consider 
a spherical nematic drop of radius a with W > 0. If 
the director is everywhere normal to the surface, as the 
surface alignment energy favors, the Frank elastic energy 
is 8ttKcl, and the surface alignment energy is zero. On 
the other hand, if the director is parallel throughout the 
interior of the drop, the Frank energy is zero, and the 
surface alignment energy is 8ttWo, 2 /3. The surface en- 
ergy scales as a 2 , whereas the elastic energy scales as a. 
Thus, surface energy dominates over elastic energy for 
large drops, and we may assume, to a good approxima- 
tion, that the preferred direction of surface alignment is 
imposed as a constraint. On the other hand for small 
droplets, elastic energy dominates over surface energy, 
and we should expect the surface director to deviate from 
its preferred orientation. The characteristic droplet di- 
mension beyond which we may assume rigid boundary 
conditions is r c = K/W « 10~ 6 /3 x 10~ 2 « 0.3^m. 
Typical droplet radii in the experiments of Poulin et al. 
p2[ are larger than 1/mi, and we may use rigid boundary 
conditions to interpret them. 

Similar considerations apply to shape distortions of the 
droplets. The positive surface tension favors spherical 
drops of either liquid crystal in water or of water in liq- 
uid crystal. The surface energy scales as aa 2 . Thus, we 
can expect drops to be spherical and undistorted by the 
nematic director for drops larger than r a = K/a « lnm. 
In what follows, we will assume that water droplets re- 
main spherical and that normal boundary conditions are 
rigidly imposed at nematic-water interfaces. 



The only topologically stable disclinations have winding 
number 1/2. All director configurations on a loop can 
either be shrunk continuously to a single point in the or- 
der parameter space, in which case the loop encloses no 
defect, or they can be continuously distorted to a path 
in RP 2 starting at some arbitrary point and ending at 
a diametrically opposite point, in which case the loop, 
encloses a disclination of strength 1/2. Typical director 
configurations for a strength 1/2 disclination are shown 
in Fig. | 





FIG. 1. Two director configurations for a strength 1/2 
disclination. In a two-dimensional nematic, the left and right 
figures correspond, respectively, to disclinations of strength 
+1/2 and —1/2. In a three-dimensional nematic, these con- 
figurations can be converted into each other via continuous 
transformations of the director 

In the simplest disclination configurations shown in 
Fig. [j], the director is n = (cos <f>/2, ± sin^/2, 0), where 
(j) = ta.n~ 1 y/x is the azimuthal angle in the aiy-plane. 
The energy per unit length of such disclination lines cal- 
culated from Eq. (0) is 



e=frK]n(R/r e ) 



(6) 



where R is the sample radius, r c is the radius of the discli- 
nation core, and e c is the core energy per unit length, 
which is of order K . 

Hedgehogs are point defects characterized by an inte- 
ger topological charge q specifying the number of times 
the unit sphere is wrapped by the director on any surface 
enclosing the defect core. An analytical expression for q 



C. Topological Defects 

Topological defects in ordered media are singular re- 
gions of spatial dimension less than that of physical 
space that are surrounded by order-parameter configu- 
rations that cannot be transformed to a homogeneous 
ground state via continuous deformations. There are two 
kinds of topological defects in a nematic. They are (1) 
line defects, called disclinations, with winding number 
of strength 1/2 in which the director undergoes a rota- 
tion of 7T in one circuit around any one-dimensional path 
encircling the linear defect core, and (2) point defects, 
called hedgehogs, in which the director sweeps out all 
directions on the unit sphere S 2 as all points on any two- 
dimensional surface enclosing the defect core are visited. 



— / dSi€ ljk n ■ (djii x d k n), 



(7) 



where the integral is over any surface enclosing the defect 
core. For an order parameter with O3 (vector) symmetry, 
the order-parameter space is S 2 , and hedgehogs can have 
positive or negative charges. Nematic inversion symme- 
try makes positive and negative charges equivalent, and 
we may, as a result, take all charges to be positive. 

There is a continuous infinity of director configurations 
for each value of the hedgehog charge. In the simplest 
unit-charge hedgehog configuration shown in Fig. ^(a), 
the director points radially outward from the point core 
like the electric field near a point charge. This config- 
uration is called a radial hedgehog for obvious reasons. 



3 



Other configurations can be obtained from the radial con- 
figuration via rotations through arbitrary angles about 
any axis. Two examples are show in Figs. ||(b) and (c). 
When the director of a radial hedgehog is rotated about 
a fixed axis through n, a hyperbolic hedgehog shown in 
Fig. ||(c) is produced. The hyperbolic hedgehog can be 
obtained from a radial hedgehog via a series of continuous 
distortions of the director passing through intermediate 
configurations such as the "circular" configuration shown 
in Fig. H (b). Thus, radial, hyperbolic, and all interme- 
diate hedgehogs are topologically equivalent. 




FIG. 2. (a) A radial hedgehog in which the director points 
radially outward from a central point like the electric field of a 
point charge, (b) A circular hedgehog obtained from a radial 
hedgehog by rotating the director at every point through n/2 
about the vertical axis, (c) A hyperbolic hedgehog obtained 
from the radial hedgehog by rotating the director at every 
point by 7r about the vertical axis. In each case, the figures 
at the left provide a three-dimensional depiction of the defect 
whereas that at the right shows a projection onto any plane 
containing the polar axis. In (b), standard notation in which 
the nail heads indicate the end of the director coming out of 
the plane is used. 

The energies of the simple hedgehog configurations 
shown in Fig. ^| in a sphere of radius R with free bound- 
ary conditions at the outer surface are easily calculated 
from the Frank free energy [Eq. (Q)]. The Frank director 



for these configurations are n = (x, y, z)/r for the radial, 
n = (y, —x, z)/r for the circular, and n = (—a;, —y, z)/r 
for the hyperbolic hedgehogs, where r = (x,y,z) and 
r = |r|. In a spherical region of radius R, their respec- 
tive energies are 

^radial = M^i - K 24 )R 

-> &tt(K - K 2i )R 

Scire = ^(3^3 + 5X 2 + 2K t - 5K 24 )R 
15 

^^L(2K-K 2i )R 

Q 

E hypcr = ^r(3^i + 2K 3 + 5K 2i )R 

^^L(K + K 24 )R, (8) 

where the final expressions are for the case of equal elas- 
tic constants. When K 24 = 0, these energies reduce 
to those calculated in Ref. Q. The hyperbolic hedge- 
hog has lower energy than the radial hedgehog provided 
K3 < 6K~i — 10/^24 or K > 2K 2 4 for single elastic 
constant approximation. Thus, if K 2 ^ = 0, the hyper- 
bolic hedgehog always has the lower energy. The circular 
hedgehog has the most bend. Since K 3 is generally the 
largest elastic constant, the circular hedgehog generally 
has the highest energy provided K 2 4 is not too large. If 
K = 2K 24 , the energies of the three hedgehog config- 
urations are equal (and equal to AttKR), as one could 
have predicted from Eq. (|2|), which is invariant with re- 
spect to rigid rotations of even a spatially varying n when 
K = 2K 2 4. In confined geometries, the Rapini-Papoular 
surface energy competes with the K 2 ^ surface term to 
determine defect configurations. 

In systems with vector symmetry, the combined topo- 
logical charge [i.e., the charge obtained by evaluating Eq. 
(0) on any surface enclosing both hedgehog cores] of two 
hedgehogs with respective charges qi and q 2 is simply 
the sum qi + q 2 . In nematics, the sign of the topologi- 
cal charge has no meaning, and the combined topological 
charge of two hedgehogs is either \q\ + q%\ or \qi~q 2 \. It 
is impossible to tell with certainty which of these possi- 
ble charges is the correct one by looking only at surfaces 
enclosing the individual hedgehogs. 

We will be primarily interested in how two unit-charge 
hedgehogs can combine to give a hedgehog charge of zero. 
Figure |3| shows how a radial and a hyperbolic hedgehog 
can combine to give a charge-zero configuration, i.e., a 
configuration in which the director is parallel at infinity. 
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radial 

hyperbolic 

FIG. 3. A radial and a hyperbolic hedgehog combining to 
give a configuration with hedgehog charge zero 

Disclination rings can carry a hedgehog charge q as 
measured by the integral in Eq. (Q) evaluated over a sur- 
face enclosing the ring |47|,|i| . Figure || depicts discli- 
nation rings with far-held director configurations corre- 
sponding to radial and hyperbolic charge 1 hedgehogs. 
These rings can be shrunk to a point leaving a point 
hedgehog. Since the disclination ring is topologically 
equivalent to a hedgehog, one can ask whether it is ener- 
getically favorable for a point hedgehog to open up to 
a disclination ring [^9|,^5). If one assumes that order 
parameter configurations remain uniaxial, one can ob- 
tain a crude estimate of the radius Rq of the disclination 
ring using the expressions, Eq. (|J) and (g), for disclina- 
tion and hedgehog energies. The director configuration 
of a charge 1 disclination ring is essentially that of simple 
disclination line discussed above Eq. (|^) in the vicinity 
of the disclination core, i.e., up to distances of order Rq 
from the ring center. Beyond this radius, the director 
configuration is approximately that of a hedgehog (ra- 
dial or hyperbolic). Thus, we can estimate the energy of 
a disclination ring of radius Rq centered in a spherical 
region of radius R to be 

B ring « 2KR Q [\nK\n{R /r c ) + e c ] + 8naK(R - R ), 

(9) 

where a = 1 — /c24 for a radial hedgehog and a = (1 + 
&24)/3 for a hyperbolic hedgehog, where fc 2 4 = K 2 a/K. 
Minimizing over Rq and setting e c = K, we find 



Ro = r c exp 



16 

7T 



7T 
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(10) 



Though admittedly crude, this approximation gives a re- 
sult that has the same form as that calculated in Refs. 
p9| , ^5|JoO| using a more sophisticated continuous ansatz. 
It has the virtue that it applies to both radial and hyper- 
bolic far-field configurations. It predicts that the hedge- 
hog with the lower energy far-field configuration (i.e., the 
one with smaller a) will have the smaller disclination-ring 
radius. If fc 2 4 = 0, the hyperbolic hedgehog has the lower 
energy with a = 1/3 rather than a = 1. In this case, the 
core of a radial hedgehog should be a ring with radius 



of the hyperbolic hedgehog, on the other hand will be a 
point rather than a ring because i?o ~ r c e~ ' 58 < r c . 

If the constraint that the tensor nematic order pa- 
rameter Qij be uniaxial is relaxed, then the core of a 
disclination can become biaxial |5l| | with a core radius 
of order the biaxial correlation length The energy 
of a disclination is still given by Eq. (|l0|) with r c ~ 
and with a core energy determined by the energy differ- 
ence between the biaxial and uniaxial state rather than 
the energy difference between the isotropic and nematic 
states. A hedgehog can also develop a biaxial core with 
radius of order £&. Because the biaxial core is charac- 
terized by a non-vanishing biaxial order parameter, its 
structure is not the same as that of the uniaxial disclina- 
tion ring discussed above. Calculations |^2| based on the 
Landau-de Gennes free energy for a nematic predict a 
biaxial core size of order 0.025/im for MBBA. A detailed 
analysis of the competition between a biaxial core and a 
biaxial disclination ring has not been done. 

It is very difficult to predict with certainty what the 
core structure of a hedgehog will be. If the core is a 
disclination ring, its radius varies exponentially with the 
elastic constants. If the core is biaxial, it will have a 
biaxial structure out to a radius of order the biaxial cor- 
relation length, which should be of order 100A or less. 
The general arguments given above would lead one to 
expect hyperbolic hedgehogs to have the smallest core 
size. In the experiments of Poulin et al. [ ^2|j34| , all hy- 
perbolic hedgehogs that were observed have cores that 
are point-like to the resolution of the optical microscope. 





Ro 



or Rq 0.2/im for 



100A The core 



(a) (b) 
FIG. 4. Disclination rings with unit hedgehog charge: (a) 
radial hedgehog and (b) hyperbolic hedgehog. The dotted line 
in each figure represents a sphere of radius R beyond which 
the director configuration is that of a hedgehog. 



III. DIRECTOR CONFIGURATIONS IN 
INVERTED NEMATIC EMULSIONS 

In the experiments reported by Poulin, et al. p2| ], a 
nematic liquid crystal (pentyl cyano biphenyl, or 5CB), 
a surfactant (sodium dodecyl sulfate, SDS) and water 
are mixed together to produce inverted and multiple liq- 
uid crystalline emulsions. The inverted emulsions are 
placed in a thin rectangular cell of approximate dimen- 
sions 20/im x 1cm x 1cm. The large-area upper and 
lower surfaces were treated to produce tangential bound- 
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ary conditions. Thus, the total hedgehog charge Q in 
the cell, obtained by performing the integral in Eq. (0) 
is zero. With normal boundary conditions, each water 
droplet nucleates a radial hedgehog of charge one. To 
maintain zero charge in the cell, compensating director 
distortions, usually point or line defects, must be created 
out of the nematic itself. Possible director configura- 
tions of a single droplet with total charge zero are shown 
in Fig. ||. A single droplet could nucleate a compan- 
ion hyperbolic hedgehog (Fig. ||a), or it could nucleate a 
disclination ring of finite radius lying above or below the 
droplet (Fig. |b) or encircling the droplet in a "Saturn- 
ring" configuration (Fig. ^|c). Director configurations for 
many droplets can be constructed from the single director 
configurations shown in Fig. |[ Other configurations in 
which the hedgehog charge of water droplets is canceled 
by continuous textures in the surrounding nematic rather 
than by the formation of point hedgehogs or singular 
disclination rings are possible. For example, if there are 
two droplets, the radial configuration around one droplet 
could continuously deform to a hyperbolic configuration 
passing through intermediate configurations such as the 
"circular" hedgehog of in Fig. ||b. The final hyperbolic 
configuration could combine with the radial configura- 
tion of the neighboring droplet to produce a configura- 
tion with zero charge but without any point defects in 
the nematic as shown in Fig. ||a. Alternatively, there 
could be a more symmetric configuration with a toroidal 
"escaped strength one" non-topological disclination line 
f33f as shown in Fig. |]b. 




(c) 



FIG. 5. Possible director configurations induced by a single 
spherical droplet with homeotropic boundary conditions in a 
nematic with total topological charge of zero, (a) Dipole con- 
figuration with a companion hyperbolic hedgehog (indicated 
by an arrow), (b) Dipole configuration with a companion hy- 
perbolic disclination ring, (c) Quadrupolar saturn ring config- 
uration with a disclination ring encircling the water droplet 
at its equator. The direction of the topological dipole p is 
shown in (a). 



radial 




hyperbolic 



circular 



(a) 




(b) 

FIG. 6. Schematic representation of the nonsingular di- 
rector configuration produced by two water droplets whose 
boundary conditions produce radial hedgehogs, (a) The ra- 
dial hedgehog around one droplet converts continuously to a 
hyperbolic configuration which then combines with the radial 
configuration of the other droplet, (b) The radial configu- 
ration around each droplet converts smoothly to a toroidal 
"escaped strength one" non-topological disclination that en- 
circles the axis defined by the droplets. We are grateful to 
R.B. Meyer for suggesting configuration (b) to us. 

In the experiments of Poulin et al. J2^,|34| , the dipole 
configuration shown in Fig. ^|a is almost always observed. 
When many droplets are in the cell, each droplet forms 
a dipole with a companion hyperbolic defect so the total 
charge of the multiple droplet system is zero as required. 
Furthermore, the droplet-dipoles align in chains parallel 
to the cell-director as shown schematically in Fig. ^ (See 
also Figs. 8 and 9 of Ref. Occasionally, droplet 

pairs are observed to induce director configurations that 
cannot be interpreted in terms of companion hedgehog 
defects [See Fig. 12 of Ref. p4[]. These configurations 
may be of the type shown in Fig. ||b. 
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FIG. 7. Schematic representation of a chain of three water 
droplets in a cell with parallel boundary conditions at infinity. 
Each droplet creates a companion hyperbolic hedgehog, and 
droplets and companions defects lie on a single line. 

In multiple emulsions, water droplets are confined to 
the interior of nematic drops with spherical symme- 
try. If the outer surface of the nematic drop enforces 
homeotropic boundary conditions, then the total topo- 
logical charge in the nematic drop is one. If there are no 
water droplets in the nematic drop, there must be a point 
hedgehog defect in the interior of the drop. In general, 
the radial hedgehog favored by homeotropic boundary 
conditions at the outer surface does not have the lowest 
energy. Instead, there is an evolution away from the ra- 
dial configuration with distance from the droplet surface 
Q,^3|, as depicted in Fig. |^a. Under crossed polarizers, 
this configuration will appear as rotating cross. A single 
water droplet in the interior of the nematic drop will cre- 
ate a radial hedgehog. Since the total topological charge 
of the nematic drop is one, no compensating defects must 
be created from the nematic. The configurations enforced 
at the water droplet surface and at the outer surface of 
the nematic drop are both radial. As a result, the director 
adopts a radial configuration throughout the drop as de- 
picted in Fig.|b (See Fig. 13 of Ref. [|§). Under crossed 
polarizers, this configuration will appear as a rigid unro- 
tated cross. A second droplet added to a nematic drop 
creates an additional interior radial hedgehog. In order 
to satisfy the global boundary condition of charge one, a 
hyperbolic defect is created out of the nematic. If there 
are N water droplets inside a nematic drop, N — 1 hy- 
perbolic defects will be created. The droplets and defects 
form linear chains with an unpaired droplet as shown in 
Fig. | (See also Fig. 14 of Ref. Q). These chains (or 
the single water droplet if that is all there is) are rigidly 
placed at the center of the nematic droplet and undergo 
no observable Brownian motion. 




(a) (b) 



FIG. 8. (a) Schematic representation of the director con- 
figuration of a nematic drop with no interior water droplet, 
the director is forced by boundary conditions to have a ra- 
dial configuration at the out surface. As distance from the 
surface increases, the director seeks lower energy, nonradial 
configurations, (b) Schematic representation of the director 
configuration of a nematic drop with a single interior water 
droplet. Homeotropic boundary conditions at the outer and 
water-droplet surfaces force a radial configuration everywhere. 




FIG. 9. A nematic drop with 3 internal water droplets. The 
3 water droplets and their 2 companion hyperbolic defects 
form a linear chain at the center of the nematic drop. The 
total charge of this configuration is one. 

IV. CONFIGURATION AND ENERGY OF 
SINGLE DROPLET 

In the preceding section, we discussed possible director 
configurations induced by the presence of spherical water 
droplets with homeotropic boundary conditions in a ne- 
matic with parallel boundary conditions at infinity. Ex- 
periments show that the water droplets create companion 
hyperbolic hedgehogs rather than disclination rings. In 
this section, our goal is to calculate the equilibrium sepa- 
ration of the droplet from its companion and to compare 
the energy of the dipole configuration with that of the 
saturn ring and intermediate configurations depicted in 
Fig. ^| The calculational program is in principle quite 
simple: solve the Euler-Lagrange equations for the di- 
rector arising from the minimization of the Frank free 
energy [Eq. (|^)] subject to the normal boundary con- 
ditions at the surface of the water droplet and paral- 
lel boundary conditions at infinity. Unfortunately, the 
Euler-Lagrange equations are highly nonlinear, and ana- 
lytical solutions cannot be found except for a few special 
geometries and boundary conditions. We can, however, 
obtain analytical solutions for the director far from the 
droplet. Using these solutions and information, evident 
from Fig. |l about the form of director configurations 
near the droplet, we can construct variational ansatzes 
for the director that obey all boundary conditions and 
that have the desired defect structure. In this section, 
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we will first discuss the nature of the far-field solutions. 
We will then use two different ansatzes to calculate di- 
rector configurations and their associated energies. The 
first ansatz applies only to the dipolar configuration. The 
second applies to all of the configurations in Fig. |jj and 
will allow us to compare, for example, the energies of the 
dipolar and Saturn-ring configurations. 

A. Far-field Solutions 

The constraint of zero topological charge requires n(r) 
to approach n = (0, 0, 1) as r — > oo. We assume that 
no is along the positive z axis. No physical results will 
change, however, if we reflect n to be along the negative 
z axis. At large but not infinite r, the deviation of n(r) 
from no is small, and n(r) « (n x , n y , 1). Thus, at large 
r, we can replace the full nonlinear Frank free energy by 
the harmonic free energy 

F hal = \K J2 J d?r{Vn^)\ (11) 

H=x,y 

where we introduced the notation n M , /i = x, y for 
the components of n perpendicular to no. The Euler- 
Lagrange equations arising from this equation are simply 
Laplace equations: 

V 2 n p = 0. (12) 

At large r the solutions to this equation can be expanded 
in multipoles: 

= — + + ^r 1 + ■■■■ ( 13 ) 

The solutions we seek are invariant with respect to ro- 
tations about the z axis and have no azimuthal compo- 
nent to n (i.e., no twist in n about the z axis). This 
implies that = and that n x and n y must be propor- 
tional, respectively, to x and y. In addition, the dipolar 
part should change sign if the position of the companion 
defect is shifted from above to below the droplet. The 
requirements are met by setting p M = (p • no)e M and 
c ij = c ( n oi e j + e i n oj) where = <5f is the unit vec- 
tor pointing in the /i = x,y direction. We identify the 
vector p as the dipole moment of the droplet-defect con- 
figuration, and p • no with its z component. Thus if p 
changes sign relative to no, the dipole contribution to n M 
also changes sign. In the configurations we consider in 
this section, p is aligned either parallel or anti-parallel 
to no so that p • no = ±p where p is the magnitude of 
the dipole moment. The parameter c, as we will show 
in more detail in the next section, is the amplitude of 
the quadrupole moment tensor cy of the droplet-defect 
combination. Thus, we have 

x zx 
n x = p z ^T + 2c— r 

n y=P*73 +2c ^5- (14) 



By dimensional analysis, p z ~ a 2 and c ~ a 3 , where a is 
the radius of the sphere. Equations |IJ produce the far- 
field configurations of Fig. ||a if we choose p z to be posi- 
tive when the companion hedgehog is below the droplet. 
Thus, we adopt the convention that the dipole moment 
of the droplet and its companion defect points from the 
companion to the droplet. 

The multipole expansion of Eqs. (|l^) and p4| ) even- 
tually breaks down because of nonlinearities neglected 
in Eq. (^lj). We can determine the leading corrections 
by including the leading anharmonic corrections to the 
harmonic free energy. Far from the defect, we can 
set n = (n x ,n y , - n 2 ± ) « (n x ,n y ,l - \n\), where 
n \ — n-^n^. The leading anharmonic correction to Fhar 
is then 

F au = $K J d 3 r(Vni) 2 , (15) 

and the Euler-Lagrange equations with this correction 
are 

V 2 n M + \n^ 2 n\ = 0. (16) 

Using this equation, one can show that if the leading 
contribution to is dipolar, then the first correction to 
rip arising from nonlinear terms is of the form r^/r 7 . In 
other words the multipole expansion of the Laplacian op- 
erator gives the correct large r behavior up to order r~ 5 . 
Thus, we could in principle develop variational approxi- 
mations in which all of the multipole moments from order 
2 to order 5 are variational parameters. We will content 
ourselves with allowing only the dipole and quadrupole 
moments to vary. 

B. The Electric-Field Ansatz 

Any ansatz for n for the dipole configuration of Fig. |E|a 
must be normal to the water droplets at r = a, tend to no 
as r — > oo, and have a hyperbolic hedgehog at some po- 
sition along the z axis outside of the water droplet. The 
familiar electrostatics problem of a charged conducting 
sphere in an external electric field can provide the basis 
for an ansatz for n that satisfies all of these conditions. 
The electric field E is normal to the conducting sphere, 
and it tends to a constant Eo = i?o e z as r — > oo. If the 
charge Q on the sphere is large enough, there is a point 
below the sphere at which the electric field vanishes. The 
normalized electric field configuration in the vicinity of 
this point is identical to that of a unit vector in the vicin- 
ity of a hyperbolic hedgehog. Thus, we have all of the 
ingredients we need for a variational ansatz. We have 
merely to choose 

n(r)=E(r)/|E(r)|. (17) 

The electric field for the above problem is 
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E 0) a , ,2 2 r a 

— — = e z + Aa^7 

En r" 3 r 



s (r 2 e 2 - 3zr), (18) 



where A 2 = Q/(Ena 2 ) is a unitless measure of the 
strength of the electric field produced by the charge Q 
compared to the fixed external field En and a is again 
the radius of the sphere. The last term in this expression 
arises from an image dipole at the center of the sphere 
that enforces the boundary condition that E be normal 
to the surface of the sphere at r = a. For A 2 > 3, we 
find precisely one zero of the electric field at r = — zne z 
outside the sphere, where zn is the appropriate solution 
to 



|z|AV + 2a 2 = 0. 



(19) 



(For A 2 = 3, the point of zero electric field just touches 
the sphere, and for A 2 < 3, a singular ring appears on 
the surface of the sphere.) zq is the distance from the 
droplet center to its companion defect. At large r, n(r) 
becomes 



(Aa 



|2 ' M 



3a J 



(20) 



in agreement with Eq. (|TJ). Thus the dipole moment 
is A 2 a 2 and the quadrupole moment is 3a 3 /2. The vari- 
able A is a variational parameter that determines both 
the position of the hyperbolic defect and the magnitude 
of the dipole moment. The ansatz fixes the quadrupole 
moment independent of the value of A and constrains the 
dipole moment to be greater than 3a 2 . A natural en- 
ergy scale is Uq = TrKa/2. The reduced energy U/Uq 
calculated from Eqs. (p"7j), (|l8|), and (|I]) is plotted as 
a function of the distance between the sphere and the 
companion defect in Fig. [H^. The energy at the mini- 
mum of the curve is U = 9.00C/o- At this minimum, the 
other parameters characterizing the droplet-defect pair 
are zq — 1.19a, p z — 3.02a 2 , and c = 3a 2 /2. 




FIG. 10. Energy (in units of Uo = irKa/2) of droplet dipole 
as a function of the distance (in units of the droplet radius a) 
from the droplet center to the companion hedgehog. 

Of course we have no a priori reason to expect a nor- 
malized electric field to yield the true minimum energy 
configuration for this problem. By modifying our cho- 
sen E(r) to be a slightly more general vector field and 
by no longer insisting that it be a true electric field, we 
may hope to improve our ansatz and to relax the con- 
straints put upon p z and c. We introduce below one such 



generalization that provides a significant improvement: 
U = 7.87U , zn = 1.26a, p z = 2.20a 2 , c = -1.09a 3 . In 
particular, we note here that the sign of c is the opposite 
of what we had previously constrained it to be through 
our electric field ansatz. 

We now present the generalization of the electric field 
ansatz E(r) in Eq. ( |l8| ) used to calculate the the re- 
sults discussed above and mention a few of the important 
properties it possesses. We insist that the generalization 
maintain the correct far-field behavior to quadrupole or- 
der so that we can identify the quadrupole strength c, as 
well as the dipole strength p z . We would also like the gen- 
eralization to put fewer constraints on the allowed values 
of p z and c. In particular, since in certain situations the 
sign of c seems to be an important quantity we should 
certainly not restrict c to be of a particular sign in the 
ansatz. Rewriting Eq. ( |l8| ) as, 



E(r) 

E 



= 1 " e, + A 2 a 2 ^ 



3za d 



(21) 



motivates the following generalization containing the 
aforementioned desired properties, 



E(r) 

En 



1 



e z + A 2 a 2 



(3ia 3 , (3 2 a k 



r k 2 +2 



(22) 

As before p z = A 2 a 2 and now c = /?ia 3 /2, provided that 
in carrying out the minimization over the appropriate 
free parameters we find k\ and k 2 are both greater than 3 
(if this had not been the case then the far-field behavior 
would not have been correct). Minimizing we find an 
energy considerably lower than that found using Eq. (fl£|). 
We obtain the results for U, zq, p z and C given above. 
In addition, we find k\ = 4.88, k 2 = 3.73, 2 = 4.27. 



C. A Second Dipole Ansatz 

To study the transition from a dipole to a Saturn ring 
and to establish that the dipole has a lower energy than 
the Saturn ring, we need an ansatz that allows disclina- 
tion rings and a limiting hyperbolic hedgehog. To con- 
struct our ansatz we use appropriately symmetric solu- 
tions to a related 2D problem and modify their far-field 
behavior to match the required 3-D far-field behavior (fol- 
lowing a route analogous to that in |35|,|36| for the case 
of the equatorial ring). 

To make contact with 2D configurations, it is conve- 
nient to look at the general problem via the following 
parameterization, 

n = (sin 0(r) cos $(r), sin O(r) sin $(r), cos O(r)) 
r = (r sin 8 cos <f>, r sin 9 sin <fi, r cos 8) (23) 

where we have expressed r in the usual spherical coordi- 
nates. The full Euler-Lagrange equations in O, $ arising 
from the free energy of Eq. (0) are simply, 
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V 2 9 -sine cos 0(V$- V*) = 

sin ©v 2 $ + 2 cos e(ve • v$) = o 

The nonlinearity of the above coupled partial differential 
equations (PDE's) make closed-form solutions difficult 
to obtain, though we note here for completeness that 
= $ = is indeed a solution (the radial hedgehog) 
and thus that solutions are not impossible to find. 

Turning now to the problem at hand, we should 
certainly impose the condition of azimuthal symmetry, 
namely: d^Q = 0, 9^$ = 1, and Q\e=o,Tr is either or ir. 
We will impose a more stringent constraint on <E>, namely 
that $ = <fi, which allows us to use our knowledge of 
2D nematic configurations to construct 3-D azimuthally 
symmetric configurations (this more stringent condition 
was also satisfied in the electric field ansatz). This leaves 
us with the single Euler-Lagrange equation, 



v 2 e- 



sin26 
2r 2 sin 29 



= 0, 



(24) 



whose nonlinearity still makes solutions difficult to ob- 
tain. Using now that we want n — > (0,0, 1) (or O — > 0) 
as r — > oo, we can linearize ( pi| ) and find the form of the 
far-field solutions, 
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At 



' J <yk~\- 

k=l 



-Pk (cost 



; Pz 



3 sin 29 



sin 9 
sin 9 sin 29 



(25) 



where the last line defines p z and c to match the defini- 
tions of these quantities in Eq. (Jl4|) , as one can check by 
substituting this form for into Eq. (p3|) . 

Given the restriction, $ = <f>, we note that in the x — z 
plane we have, 



n2D = (sin 0, cos 0). 



(26) 



For ri2_D, we have a wealth of information about how to 
construct solutions for the harmonic free energy, 



Fid = \K 



d 2 r(V 2 D@) 2 - 



(27) 



It would be nice if our problem reduced to this linear 
problem. However, even though this is not the case we 
can still exploit our knowledge of the 2D solutions to 
construct ansatz solutions for the 3D problem. The pro- 
cedure is quite simple and has at least some promising 
motivations. Any configuration of ii2D that is invariant 
under x — > —x can be converted to a 31? configuration 
by spinning about the z axis to produce 



Furthermore we note that the canonical q 2 D = +1 de- 
fect on the z— axis becomes a q = 1 radial hedgehog in 
3D, and the canonical q 2 D = —1 defect on the z— axis 
becomes a q = 1 hyperbolic hedgehog in 3-D (recall that 
in 2D the charges of nematic defects are signed and their 
composition law is addition). Also, any symmetric pair 
(we need a pair to maintain the required reflection sym- 
metry about the z— axis) of q 2 o — if defects off the 
z— axis, when spun into a 3-D configuration become a 
strength i— disclination ring. Finally, we note that sat- 
isfying the 3-D BC on the sphere, namely that n = e r on 
the sphere, simply requires satisfying normal BC on the 
circle in 2D. In 2-D we can satisfy these BC using the 
method of images, which works because of the linearity 
of the EL equations. 

Before writing down the ansatz we note that the 2-D 
configuration = 9 ZQ , where 9 Za is the polar angle 
measured with respect to the point zq on the z— axis, 
when made into a 3-D configuration indeed satisfies the 
3D Euler-Lagrange (EL) equations. In two dimensions, 
©i + 02 is a solution to the EL equations provided both 
©i and ©2 are. However, because the 3D EL equations 
are nonlinear, ©i + ©2 is in general not a solution to the 
3-D El equations even if ©i and ©2 individually are. It 
is this fact that prevents our ansatz configurations from 
being true solutions to the 3-D equations. 

We now construct our ansatz for a sphere at the origin 
with a compensating unit strength hyperbolic hedgehog. 
To keep equations as simple as possible, we use units in 
which the sphere radius is one. As discussed above, we 
first construct a solution to the two-dimensional problem. 
This is done using the fact that a defect of strength q at 
position r 2 _D = (xo,zo) is described by the field = 
q , tan _1 [(2; — zq)/(x — xq)]. Boundary conditions at the 
sphere's surface and at infinity can be met by placing 
a strength q = +2 defect at the sphere's center, a — 1 
defect at position (0, — r^) for arbitrary r^, and an image 
— 1 defect inside the spheres at (O,^ 1 ) as shown in Fig. 
[Tlla,. This leads to 



©o = 29 - tan" 1 



r sin 9 



tan 



rrd sin t 



rcost 



rrd cos 9 + 1' 



(29) 



where we have taken the radius of the sphere to be 1, and 
ra is the distance from the defect below the sphere to the 
origin. While this does have the correct form near the 
defect at z = — rd it does not have the correct far field 
form for the 3D problem. In fact, 



n = (sin © cos (f>, sin sin cf>, cos 0) , 
a solution with the $ = <b constraint. 



(28) 
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(a) (b) (c) 

FIG. 11. Utilizing the method of images we are able to construct solutions for the related 2D problem of a nematic with 
homeotropic boundary conditions at infinity and with normal boundary conditions on a circle about the origin. These config- 
urations are then extended to 3D configurations by spinning them about their vertical symmetry axis, where the singularity 
along the vertical axis become singularities in 3D (possibly removable ones as in the case of the +2 defect at the origin) and the 
symmetric pairs of defects off the axis become singular defect rings. Note that in our calculations we are not really concerned 
with the from of the field inside the circle since the nematic is only present in the exterior, the field is merely drawn here to 
elucidate the origin of the ansatz used, (a) a hyperbolic defect beneath the sphere, (b) a non-equatorial disclination ring, (c) 
an equatorial disclination ring. 



e 



i 

r,i 



1 \ sin 9 cos 9 



(30) 



which does not agree with Eq. (|25|). However, we note 
that to the order shown, it differs only by an overall power 
of r _1 (this is not true for the higher order terms not 
shown). So, we alter 0o in the following way, taking care 
to preserve the BC at r = 1. 



and U = 7.84J7o in very good agreement with the results 
obtained from the generalized electric field ansatz. In 
addition r<2 = 1.22a agrees nicely with the corresponding 
quantity zq — 1.26a obtained from the electric electric- 
field ansatz. Preliminary results obtained via numeri- 
cal solution to the full Euler-Lagrange equations for this 
problem arc in excellent agreement with these ansatz re- 
sults 



6= 29 - tan" 1 

-k/r 



r cos 9 + rd 
1 



tan 



1 1 
- + ^ 



rr<j sin t 



rr d cos 9 + 1 



4) {h-h Uwl ° 



which now has the far-field form, 
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r d + — 



smt 



sin26> 



I'd 



(31) 



(32) 



and 



This agrees with Eq. ( |25| ) with p z 
c = — (r 2 , + r~ d 2 )/2, which contrary to the electric field 
ansatz is negative (recall a = 1). The factor e~ fe / r3 , in- 
troduced for numerical convenience, tends to 1 at large 
r and to a value near the sphere controlled by the varia- 
tional parameter k. Substituting this form [Eq. (|3l])] for 
(and $ = 0) in Eq. (^|) and minimizing over Td and k 
(numerically) we find, k — 0.32 and — 1.22. 

Now taking the sphere to have a radius of a, as in the 
electric-field ansatz, we find p z = 2.04a 2 , c = —1.08a 3 



D. From Dipole to Saturn Ring 

We can easily generalize the angular parametrization 
of a point defect just discussed to describe an annular 
ring defect with a varying opening angle 9d [Fig. |i~l|b] and 
thereby study the transition from a dipole configuration 
wth 9d = 7T to the saturn ring with 9d — n/2 [Fig. [Tl|c]. 
We proceed exactly as in the point-defect case. We place 
one strength +2 defect at the center of the circle, two 
strength —1/2 defects at r2D — sm9d, cos 9d), and 
two strength —1/2 images inside the sphere at Y2D = 
r~ l 1 (±sm9d,cos9d). This gives a solution to the 2D 
problem with two strength —1/2 defects outside the cir- 
cle. When promoted to 3D, this solution correctly sat- 
isfies homeotropic boundary conditions at the surface of 
the sphere, and it yields a strength 1/2 disclination ring 
outside the sphere with an opening angle of 9d- It fails, 
however to produce the correct far-field form for n. We 
add terms similar to those of the preceding calculation 
to correct this deficiency to produce 



0=29-- 
2 



tan 



rsmf 



Td sin ( 



- tan 



r cos 9 — Td cos 9d 
r sin 9 + rd sin^d 
r cos 9 — rd cos 9 d 
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_ i rrd sin 9 — sin 64 _ 1 rrd sin 9 + sin 9d 

tan t — h tan 



-fc/r 3 



rr^ cos t/ — cos Ud 

1 \ /l 1 



rr^ cos 6> — cos &d 



cos t/rf sin ( 



^+4) (^-^) ( "" i - 2< '"^ w ' ,!sii,w<usf 



1 \ / 1 



cos 6^(4 cos 2 6>d - 3) x 



sin (9(4 cos 2 6» - 1) 



(33) 



The director configuration in the vicinity of the discli- 
nation ring is singular, and the above form breaks down 
at distances from the disclination ring less than the core 
radius r c . The equatorial ring configuration, 9d = tt/2 
has been investigated previously |55|,[56| . 

Figure |lj shows the energy U in units of Uq = nKa/2 
for various values of 9d obtained by minimizing the free 
energy [Eq. (||)] , augmented by an additional core energy 
irrdK sin #(2/2, over the variational parameters and k 
in the ansatz function Eq. (33). The reduced core ra- 
dius in these calculations was chosen to be r c = 10~ 3 . 
We have checked that our results are insensitive to the 
value of r c for 10 -4 < r c < 1CT 2 . For a typical core ra- 
dius of lOnm, or results are thus good for particles with 
radii as small as 1/im. Figure [l^ shows the correspond- 
ing equilibrium distance of the ring from the center of the 
sphere. We note that Td and E for 9d ~ 7r agree quite 
well with the values obtained from the point defect ansatz 
previously presented. That is, this ansatz does collapse 
down to the point defect in a nice manner. However, we 
must also note that for 9d — the equatorial ring, our 
rd = 1.08a is somewhat different from the value of 1.25a 



found by Terentjev in 35 3 




FIG. 12. Energy versus angular position of the defect ring. 
Note that the equatorial ring (0<j = j) does appear to enjoy 
some metastability but that the collapsed ring (or effectively 
the point defect, — n) is of much lower energy. 




nil 3;t/4 7t 

FIG. 13. The preferred distance from the origin of the 
disclination ring as a function of its angular position. 

Further comments about the ring configuration are in 
order. Recalling that in three dimensions the point defect 
singularities represent integrable singularities whereas 
the ring defects do not, one might naively expect that 
ring s sho uld always collapse into points. As discussed in 
Sec. IIC, however, this is not always the case. Even an 



isolated point singularity might have higher energy than 
an isolated ring singularity |^,^5|,^0| . But it is worth 
noting that in such circumstances the equilibrium radius 
of such rings turns out to be quite small (« 0.2/im). Thus 
here we would naively expect that if a ring configuration 
were to exist it would probably not be an equatorial ring 
configuration. This proves to be correct as we see in 
Fig. |l2]. Though the equatorial ring does appear to en- 
joy some metastability, its energy is considerably higher 
than that of the point defect below the sphere. 



E. Thermal Stability 

We have seen so far that the point defect beneath the 
sphere is the energetically favorable configuration. The 
elastic constant k for deviations Az from equilibrium sep- 
aration zq is simply the curvature of the energy versus 
separation curve [e.g., Fig. [l(J. The dipole ansatz yields 
k = 33nK/a. The other ansatzes yield similar values. 
Thus 




k R T 



10" 



(34) 



where the final numerical estimate follows from ksT rj 
10~ 13 erg, K « 10~ 6 dyne, and a = 1/im. These fluctu- 
ations in the length of the topological dipole are unob- 
servably small. 

We have argued that the topological dipole prefers to 
align parallel or anti-parallel to the director at infinity. 
We will now show that the angular restoring force con- 
stant kg is greater that 2irKa so that 



((A9f) < 



k B T 
4-KKa 



10" 



(35) 
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Thus, though of order ten times greater than fluctuations 
in the length of the dipole, angular fluctuations are still 
unobservably small. Interestingly, we note that angular 
fluctuations in the 2D version of this problem are much 
larger and have indeed been observed in free standing 
smectic films ]3!| 

Our approach is to provide an ansatz for a director con- 
figuration with the dipole rotated through an angle AO 
relative to the director field at infinity given a configura- 
tion in which the dipole moment is parallel to the director 
field at infinity. We will then use this ansatz to calculate 
bounds on kg. We start with an aligned dipole config- 
uration with n expressed in polar coordinates according 
to Eq. (p3|). We then construct a rotated configuration, 
n', by slowly rotating n about the y-axis as we progress 
radially outward: 



n^, = sin 6 cos <f> cos / + cos O sin / 



n' = sin sin <!> 



n' = cos cos / — sin cos $ sin / 



(36) 



where the amount we rotate at each point is given by 
the function /(r) which must satisfy the boundary con- 
ditions, 



f(r = a) = 
f(r = R) = AO. 



(37) 



AO denotes the tilt angle of the dipole with respect to the 
far-field, a the droplet radius, and R the system size. To 
see that AO is the stated tilt angle we note that the far- 
field of n' makes an angle AO with the z— axis, and fur- 
thermore, the transformation of n to n' does not change 
the position of the singularity, or of the droplet itself. 
Thus, the droplet-defect dipole (p) is still aligned with 
the z— axis. 

We denote F(A0) as the free energy of the tilted dipole 
configuration and, accordingly refer to F(0) as the free 
energy of the aligned configuration. Using Eq.(H) we find, 



F(A0) - F(0) = 



\K J d 3 r(cos 2 $ sin 2 + cos 2 0) V/ • V/ 



(38) 



where we have eliminated terms linear in / using the 
symmetry F(A0) = F(—A0). Now noting that the par- 
enthetical term is always less than unity we have, 



F(A0) - F(0) > \K J d 3 rVf • V/. 



(39) 



The right hand side of this equation is a minimum when 
V 2 / = subject to the boundary conditions of Eq. (37), 
i.e., when 



f = A0 



1 - (o/r) 



1 - (a/R)' 
Using this / in Eq. (|39]), we obtain 

F(A0) - F(0) > 2ttK{A0) 2 cl, 
for R 3> a, implying kg > AitKa. 



(40) 



(41) 



F. Optical Images 

In the experiments reported by Poulin et al. p2| , only 
the dipole and not the saturn ring configuration shown 
in Fig. H is observed. Figure [lip , presents an experi- 
mentally obtained image of a single water droplet under 
crossed polarizers with one polarizer parallel to the dipole 
axis. In the region of the droplet we see a pronounced 
pattern arising from the spatially varying director field. 
In Fig. 0b we show an image of a similar single droplet 
calculated using the Jones matrix formalism pq j and ne- 
glecting any refraction at the droplet boundary. The sim- 
ilarity of the two images is obvious and clearly confirms 
the occurrence of the dipole configuration. Both pictures 
show two bright wings left of the droplet. In the cal- 
culated picture they are much more extended than they 
are in the experimental picture. The theoretical picture 
was calculated using the simple electric field ansatz. The 
more sophisticated ansatzes reduce the region around the 
defect where there is rapid variation of the director. They 
would yield images in closer agreement with the experi- 
mentally observed one. 

FIG. 14. (a) Image of a single droplet with its companion 
defect as observed under crossed polarizers obtained by P. 
Poulin. (b) Simulated image of the same configuration using 
the Jones matrix formalism. The two images are very similar. 



V. PHENOMENOLOGICAL THEORY AND 
DROPLET-DROPLET INTERACTIONS 

To understand the properties of multi-droplet emul- 
sions, we need to determine the nature of droplet-droplet 
interactions. These interactions are mediated by the ne- 
matic in which they are embedded and are in general 
quite complicated. Because interactions are determined 
by distortions of the director field, there are multi-body 
as well as two-body interactions. We will content our- 
selves with calculations of some properties of the effec- 
tive two-body interaction. To calculate the position- 
dependent interaction potential between two droplets, we 
should solve the Euler-Lagrange equations, as a function 
of droplet separation, subject to the boundary condition 
that the director be normal to each droplet. Solving com- 
pletely these nonlinear equations in the presence of two 
droplets is even more complicated than solving them with 
one droplet, and again we must resort to approximations. 
Fortunately, interactions at large separations are deter- 
mined entirely by the far-field distortions and the multi- 
pole moments of the individual droplet-defect pairs, and 
they can be described by a phenomenological free energy, 
which we will derive in this section. 

In the preceding section, we established that each wa- 
ter droplet creates a hyperbolic hedgehog to which it 
binds tightly to create a stable topological dipole. The 
original droplet is described by three translational de- 
grees of freedom. It draws out of the nematic a hedge- 
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hog, which itself has three translational degrees of free- 
dom. The two combine to produce a dipole with six 
degrees of freedom, which can be parametrized by three 
variables specifying the position of the water droplet, two 
angles specifying the orientation of the dipole, and one 
variable specifying the magnitude of the dipole. As we 
have seen, the magnitude of the dipole does not fluctu- 
ate much and can be regarded as a constant. The di- 
rection of the dipole is also fairly strongly constrained. 
It can, however, deviate from the direction of local pre- 
ferred orientation (parallel to a local director to be de- 
fined in more detail below) when there are many droplets 
present. The droplet-defect pair is in addition charac- 
terized by its higher multipole moments. The direction 
of the principal axes of these moments is specified by 
the direction of the dipole as long as director configura- 
tions remain uniaxial. The magnitudes of the uniaxial 
moments like the magnitude of the dipole moment are 
energetically fixed. When director configurations are not 
uniaxial, multipole tensors will develop additional com- 
ponents, which we will not consider here. We can thus 
parametrize droplet dipoles by their position and orien- 
tation and a set of multipole moments, which we regard 
as fixed. Let e a be the unit vector specifying the di- 
rection of the dipole moment associated with droplet a. 
Its dipole and quadrupole moments are then p Q = pe a 
and cfj = c(e"e" — g<%), where p and c are the magni- 
tudes of the dipole and quadrupole moments calculated 
in the preceding section. We can now introduce dipole- 
and quadrupole-moment densities, P(r) and Cy(r) in the 
usual way. Let r Q denote the position of droplet a, then 

P(r) = ^p Q (5(r-r") 

a 

CyW^c^r-r"). (42) 

a 

We now construct an effective free energy for direc- 
tor and droplets valid at length scales large compared 
to droplet dimensions. At these length scales, we can 
regard the droplets as point objects (as implied by the 
definitions of the densities given above). At each point 
in space, there is a local director n(r) along which the 
droplet-dipoles wish to align. In the more microscopic 
picture, of course, the direction of this local director cor- 
responds to the far-field director no. The effective free 
energy is constructed from rotationally invariant combi- 
nations of Pi , Cij , m , and the gradient operator Vt that 
are also even under n — * — n. It can be expressed as a 
sum of terms 

F = F n +F p +F C + Malign , (43) 

where F n is the Frank free energy, F p describes inter- 
actions between P and n, Fc describes interactions be- 
tween Cij and n involving gradient operators, and 

FeXi ga = -D / dV<7y(r)n;(r)n 3 (r) 



= -#Q^[(e«.n(r«)) 2 -I] (44) 

a 

describes the alignment of the axes e a along the local di- 
rector n(r a ). The leading contribution to F v is identical 
to that for electric dipoles in a nematic j2^j5(| 

F p = AttK J d 3 [-P • n(V • n) + (3P • (n x V x n)], 

(45) 

where f3 is a material-dependent unitless parameter. The 
leading contribution to Fq is 

F c = 4ttK J d 3 r[{V ■ n)n • V(n i Cy 7i J ) 

+V(mCijTij) • (n x V x n)]. (46) 

There should also be terms in Fc like CyV 'kUiV \rij. 
These terms can be shown to make contributions to the 
effective droplet-droplet interaction that are higher order 
in separation than those arising from Eq. (|46|). Equa- 
tion (|46l ) is identical to that introduced in Ref. J37[| to 
discuss interactions between saturn-ring droplets, pro- 
vided fiiCijiij is replaced by a scalar density p(r) — 
^2 a 5(r — r Q ). The two energies are absolutely equivalent 
to leading order in the components n M of n perpendicu- 
lar to n provided all e° are restricted to be parallel to 
no. When this restriction on e a is lifted, and to higher 
order in n M , the two theories differ. In our opinion, the 
scalar variable cannot strictly speaking be used, because 
each droplet carries with it an anisotropic director envi- 
ronment, even when the dipole moment is zero. 

Since P prefers to align along the local director n, the 
dipole-bend coupling term in Eq. ^| can be neglected to 
leading order in deviations of the director from unifor- 
mity. The — P • n(V • n) term in Eq. ( p5| ) shows that 
dipoles aligned along n create local splay as is evident 
from the dipole configuration depicted in Fig. ||a. In addi- 
tion, this term says that dipoles can lower their energy by 
migrating to regions of maximum splay while remaining 
aligned with the local director. Experiments |2^,|34| sup- 
port this conclusion. Water droplets in a nematic drop 
with homeotropic boundary conditions at its surface con- 
gregate at the center of the nematic drop where the splay 
is a maximum. Boundary conditions at the outer surface 
of nematic drops can be changed from homeotropic to 
tangential by adding a small amount of glycerol to the 
continuous water phase. In the passage from homeotropic 
to tangential boundary conditions, the topological charge 
of the nematic drop changes from one to zero, and point 
defects called boojums |3(^,[25],[3lJ form on the drops sur- 
face. The director splay is a maximum in the drop's 
interior near the boojums. Water droplets move from 
the drop centers to drop surfaces near boojums as the 
boundary conditions are changed. The final configura- 
tion of two droplets in a nematic drop with tangential 
boundary conditions is shown in Fig. |l5|. 
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FIG. 15. Schematic representation of two water droplets 
with homeotropic boundary conditions at their outer surfaces 
in a nematic drop with tangential boundary conditions at its 
outer surface. The total hedgehog charge is zero, so there is 
one hyperbolic hedgehog per droplet. The water droplets mi- 
grate to the region of maximum splay near a surface boojum. 
the splay near the boojum is assumed to be sufficiently strong 
that both dipoles prefer to be near the boojum rather than 
to form a chain. 



To harmonic order in n^, the full effective free energy 



F = K J d 3 r[i(Vn M ) 2 - ^P z d^ + 4^(c9 z C zz )d p 



(47) 



The dipole-bend coupling term of Eq. ( J45| ) does not con- 
tribute because P is aligned along the far-field director. 
Thus, 



VV = H[P z (r)-a z C zz (r)], 



(48) 



or 



?v( r ) 



-X[P z (v')-d' z C zz (v')}. (49) 



p z S(r) (p z = ±p), and C Z2 (r]_= |c<5(r), and the above 



For a single droplet at the origin with e = n , P z (r) 

w = " 

equation yields exactly Eq. (|lj). 

Droplets create far-field distortions of the director, 
which to leading order at large distances are determined 
by Eq. (fl8|), that interact with the director fields of other 
droplets. This leads to an effective droplet-droplet inter- 
action that can be expressed to leading order as pair- 
wise interactions between dipole and quadrupole densi- 
ties. Using Eq. @ in Eq. @, we obtain 

^ = \j d\d?r'[P z {v)V PP {v r')P z (r>) 

+C zz {t)V C c{v-v')C zz {v') (50) 
+V PC (r - v')[C zz {v)P z {v') - P z {v)C zz {v')] 



Vpp(r)=d M c^ = ^(l-3cos 3 



1 



Vcc(r) = -d 2 d^ z = — (9 - 90 cos 2 9 + 105 cos 4 1 



Vpc(r) = dzdf.dfj.- = 



r r 
1 cos 9 



(15 cos 2 9 - 9), 



(51) 



where 9 is the angle the separation vector r makes with 
no . The interaction energy between droplets at positions 
r and r' with respective dipole and quadrupole moments 

Pz, p' z , c 



U(R) = 4ttK PzP ' z V pp (R) + -ccV cc (R) 



-(cp' z -c'p z )Vp C (R) 



(52) 



This potential can be used to calculate the force between 
two droplets as a function of their separation. Consider, 
for example, the interaction between two droplets labeled 
1 and 2 with respective radii a\ and For simplicity, as- 
sume the dipoles associated with each droplet are aligned 
along the positive z axis and that the center of droplet 
1 is at the origin and that of droplet 2 at r = (0, 0, R) 
a distance R away along the positive z axis as shown in 
Fig. The dipole and quadrupole moments scale re- 
spectively as a 2 and a 3 , and we can write p z = aa 2 and 
c = — 3/3a 3 /2. The dipole ansatz solution of Sec. IVC, 
predicts a = 2.04, nd (3 = (2/3) x 1.08 = 0.72. The force 
between two droplets is then 



ti 



AttK 



12 i?4 

,22/ 



1 20 



-aj3a 1 a 2 {a\ — 02) 



24 

w 



(53) 



The dominant force is the attractive dipole-dipole force 
proportional to i?~ 4 . Recent experiments confirm this 
relation ^7j. Interestingly the sign of the dipolc- 
quadrupole force, which dies off as i?~ 5 , vanishes for 
particles of equal radius. When the particle have un- 
equal radii, the sign of this force depends on the relative 
position of the large and small particle. If a\ < <Z2j it is 
repulsive (for > 0); if a\ > a^, it is attractive, i.e., it is 
repulsive if the smaller ball is to the right (positive z) of 
the large ball and attractive if it is to the left (negative 
z). 



where 
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FIG. 16. (a) configuration in which the large particle is to 
the right of the small particle, (b) the inverse configuration. 
The force between the two particles is more attractive in case 
(a) than in case (b). In both cases, the particle to the left is 
labeled 1 and that to the right 2 

VI. SUMMARY AND CONCLUSIONS 

Inverse nematic emulsions in which surfactant coated 
water droplets are dispersed in a nematic host have prop- 
erties that are distinct from those found in colloids, emul- 
sions of two isotropic fluids, and of emulsions of nematic 
droplets in an isotropic fluid. The water droplets in these 
emulsions exhibit anisotropic interactions that are repul- 
sive at short range and attractive at long range. The 
short-range repulsive interaction prevents coalescence of 
droplets and leads to long-term stability, which can be 
eliminated by heating into the isotropic phase. The long- 
range attractive force is dipolar and favors chaining of 
droplets. 

In this paper, we have presented a detailed theoreti- 
cal study of droplets and droplet interactions in inverse 
nematic emulsions. Homeotropic boundary conditions 
at droplet surfaces produce a hedgehog director config- 
uration around each droplet. Constraints on the global 
topological charge force the nucleation of compensating 
topological defects out of the nematic host. The com- 
pensating defect associated with a single droplet in a cell 
with a parallel aligned director at infinity can be a point 
hedgehog or a disclination ring sitting above or below the 
droplet or encircling its equator in the saturn ring con- 
figuration as shown in Fig. |^. Using various variational 
ansatzes, we showed that in the lowest energy configura- 
tion, a single water droplet pulls a single point hedgehog 
from the nematic to form a tightly bound dipole. Then, 
using a phcnomcnological model in which the topological 
dipoles are coupled to the nematic director via a flexo- 
electric interaction, we derived the effective long-range 
dipolar interaction between water droplets. We also con- 
sidered quadrupolar corrections to the dominant dipolar 
interaction. The phcnomcnological model also predicts 
the experimentally observed tendency of dipoles to seek 
regions of high splay. 

We have focussed mostly on interactions between 



droplets in cells with parallel boundary conditions at in- 
finity with total topological charge zero. Multiple emul- 
sions in which water droplets are dispersed in nematic 
drops, which are in turn dispersed in water, have made 
possible the isolation of a finite number of droplets and 
facilitated a number of experimental observations. The 
nematic drops are characterized by a topological charge 
of one rather than zero and by spatially nonuniform di- 
rector configurations. Many of the properties of these 
droplets-within-drops systems such as chaining and the 
tendency of the water droplets to concentrate near the 
center of the nematic drop are explained by the analy- 
ses in this paper. Numerically accurate predictions about 
these systems, however, require, global minimization pro- 
cedures that can only be done numerically. Numerical 
algorithms to study droplets dispersed in confined ge- 
ometries are currently under development f55fl . 

Inverse nematic emulsions are a relatively new addi- 
tion to the ever growing list of interesting soft materials, 
and they offer the hope of new and surprising proper- 
ties. We are currently investigating among other things 
the dynamics of droplets in inverse emulsions and inverse 
emulsions of water droplets in cholesteric rather than ne- 
matic liquid crystals. 
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